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Accepting that galactic mass discrepancies are due to modified dynamics, I show why it is specif¬ 
ically the MOND paradigm that is pointed to cogently. MOND is thus discussed here as a special 
case of a larger class of modified dynamics theories whereby galactic systems with large mass dis¬ 
crepancies are described by scale-invariant dynamics. This is a novel presentation that uses more 
recent, after-the-fact insights and data (largely predicted beforehand by MOND). Starting from a 
purist set of tenets, I follow the path that leads specifically to the MOND basic tenets. The main 
signposts are: (i) Space-time scale invariance underlies the dynamics of systems with large mass 
discrepancies, (ii) In these dynamics, G must be replaced by a single “scale-invariant” gravitational 
constant, Qo (in MOND, Qo = = Gao, where Oq is MOND’s acceleration constant), (in) Uni¬ 

versality of free fall points to the constant go = Qo/G as the boundary between the G-controlled, 
standard dynamics, and the Qo-controlled, scale-invariant dynamics (in MOND, go = uo)- (iv) Data 
clinches the case for go being an acceleration (MOND). 

PACS numbers: 


I. INTRODUCTION 


MOND [1, 2] is a paradigm of dynamics that attributes 
the mass discrepancies in galactic systems not to dark 
matter but to a major departure from Newtonian dy¬ 
namics and general relativity at low accelerations. The 
basic tenets of the MOND paradigm pertain to the weak- 
field limit (in the relativistic sense; potential <C c), and 
are: (i) The modified dynamics involve a new constant 
with the dimensions of acceleration, Oq. (h) Dynamics 
go to standard dynamics when system accelerations are 
much larger than Oq (formally, in the limit Oq —>-0). (iii) 
Dynamics become scale invariant (SI) for accelerations 
much below Oq, which characterize the dynamics of sys¬ 
tems that exhibit large mass discrepancies. There, New¬ 
ton’s constant, G, which obstructs SI, has to be replaced 
by the “SI” constant, Aq = OpG, as the “gravitational 
constant”. Formally, this limit is gotten by Oq —> oo, 
G —>■ 0, with ^0 fixed. References [3, 4] are recent re¬ 
views of MOND. 


My purpose here is to show, systematically, how one 
is lead to these basic tenets, almost inescapably, by the 
data on galactic mass discrepancies, once one attributes 
these discrepancies to modified dynamics. 


In Sec. II, I explain why we should accept that galac¬ 
tic systems with large mass discrepancies are governed 
by SI dynamics. Some general consequences of SI in sys¬ 
tems governed by gravity are discuss in Sec. III. Section 
IV discusses the transition from the modified regime to 
the standard regime. In Sec. V, I explain why of all 
the possible SI dynamics, it is MOND-with the pivotal 
role of accelerations-that is pinpointed. Section VI is a 
summary and discussion. 


II. SPACE-TIME SCALE INVARIANCE 

UNDERLIES LARGE-MASS-DISCREPANCY 
DYNAMICS 

We start with the observation that the circular rota¬ 
tional speeds around galaxies become asymptotically flat. 
This is evinced, in the first instance, by the measured ex¬ 
tended rotation curves of discs galaxies (see, e.g.. Figs. 
15, 21-27 of Ref. [3], and the figures in Ref. [5]). The 
asymptotic-flatness region occurs at large enough radii 
where the baryonic mass may already be taken as a point 
mass, and where we thus expect a Keplerian fall-off of the 
rotational speeds. 

Asymptotic flatness of rotation curves-with its clear 
conflict with the prediction of standard dynamics with¬ 
out dark matter (Kepler’s 3rd law)-is, arguably, the most 
clear-cut and most iconic evidence for a mass discrepancy 
in galactic systems. These asymptotic regions (and oth¬ 
ers) are thus characterized by large mass discrepancies. 

In advancing MOND [1], I have elevated this obser¬ 
vation to a cornerstone and an axiom in constructing 
modified-dynamics explanations of the mass discrepan¬ 
cies: In the thought-for modified-dynamics theory the 
circular orbital speed around an isolated, bounded mass 
becomes asymptotically constant, V{r) —>■ 

Rotation curves pertain to circular motions of massive 
test particle. Asymptotic flatness tells us that circular- 
orbit periods, around a given mass, become proportional 
to the orbital radius at radii where the mass discrepancy 
is large. In other words, at such radii all circular orbits 
around a given mass are related to each other by a scaling 
transformation 

(r,t) ^ A(r,t). (1) 

In yet other words, the family of asymptotic, circular 
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orbits around a central body is space-time Sid 

There is also evidence that such a space-time scaling 
applies to other types of orbits in the regions of large 
discrepancies, (i) The data and analysis of Refs. [6, 7] 
and the subsequent analysis of Ref. [8]-which concerns 
extended hot gas in hydrostatic equilibrium around two 
isolated, elliptical galaxies-is tantamount to showing this 
SI for the isotropic distribution of velocities in the gas. 
(ii) The weak-lensing analysis of Ref. [9] (see also com¬ 
plementary analysis in Ref. [10]) extends the evidence for 
SI to unbound photon trajectories. The analysis of Ref. 
[9], while consistent with other behaviors, is also consis¬ 
tent with the gravitational-light-bending angle becoming 
asymptotically independent of the impact parameter (for 
radii where the central galaxy is dominant).^ This im¬ 
plies that asymptotic null-geodesics can be gotten from 
each other by space-time scaling. 

All these would follow if the modified dynamics that 
describe phenomena exhibiting large mass discrepancies 
are SI. By this one means that if ri{t) are the solution 
trajectories, for a self gravitating system of masses nii, 
for some initial conditions, then, Xri(t/X) is the solution 
for the appropriately scaled initial conditions. 

A poor-man’s general statement of scale invariance, in 
the present context, is that velocities-which are invari¬ 
ant to scaling-should depend only on masses, which are 
invariant, but not on sizes, which are not. This is mani¬ 
fested in the asymptotic flatness of rotation curves, and 
would also underlie the strong mass-velocity correlations 
that are known to hold in galactic systems, such as the 
Tully-Fisher and the Faber-Jackson correlations. 

Indeed, many years after the advent of MOND it was 
realized [11, 12] that this symmetry can be taken as a 
basic tenets for the modified dynamics, and many of 
the observations that concern phenomena in regions of 
large mass discrepancies automatically follow (without 
dark matter). But MOND posits more than SI in large- 
discrepancy systems; it also pinpoints the special role of 
accelerations. It is then left to be shown how this is 
forced on us. 

I thus take here the SI in the limiting theory describ¬ 
ing high-discrepancy systems as the starting point for 
constructing the modified dynamics.^ Call this limiting 
theory the DML theory.^ 

As always in such generalizations, the principle is much 
more general than the observations on which it is based. 


^ It is assumed that the asymptotic dynamics depends only on the 
total mass of the central object and not on its size attributes, 
which also change under scaling. 

^ This is equivalent to the finding of the above reference that the 
observed light bending is consistent with a logarithmic potential. 

® By SI I mean only that the equations of motion are SI, not nec¬ 
essarily the action. For that it is enough that the action has a 
well-defined scaling dimensions, namely, that it is multiplied by 
some constant (a power of A) under scaling. 

Since in the end this theory will be identified with the deep- 
MOND-limit. 


III. CONSEQUENCES OF SCALE INVARIANCE 

We do not know whether the modification of dynam¬ 
ics invoked to explain the discrepancies in the gravity- 
governed galactic systems has bearings also on non- 
gravitational physics, such as electromagnetism. Be the 
case as may be, here I confine myself to the purely grav¬ 
itational sector of such a modified dynamics. Further¬ 
more, I confine myself to the nonrelativistic (NR) limit. 
In addition, and in harmony with this, I assume that 
in the modified dynamics, bodies (such as galaxies) can 
still be assigned the attribute of a mass, which is the 
only source of their gravity, and which coincides with the 
standard notion of mass.® This mass may have to be 
measured by means outside the DML theory, for exam¬ 
ple, its value may be assigned relative to atomic masses 
or to the Planck mass. Or, it may be taken as the sum of 
the masses of its constituents (atoms, stars, etc.), which 
can be determined with standard Newtonian gravity. ® 

A SI theory can be recognized by its constants; i.e., by 
whether they retain their values under a simultaneous 
change of the length and time units by the same factor. 
More precisely, the degrees of freedom of a theory can be 
normalized to have units such that this becomes a nec¬ 
essary and sufficient condition for SI.^ In other words, 
with this standardized choice of units for the degrees of 
freedom, which I assume henceforth with no loss of gen¬ 
erality, the constants must have dimensions of the form 

For example, the fact that Newtonian dynamics (and 


® In principle, the modified dynamics could be such that gravity 
depends on attributes of the sources other than mass. 

® When we say in the context of the mass-discrepancy problem 
that the baryonic mass falls very short of the measured dynam¬ 
ical mass, we mean by “baryonic mass” the sum of the masses 
of the system’s constituents, measured individually by various 
means. For example, we measure the baryonic mass in neutral 
hydrogen, by measuring the 21 -cm luminosity, which is thought 
to be proportional to the number of hydrogen atoms, which we 
multiply by the laboratory value of the hydrogen mass. The 
baryonic mass in stars is based on masses of individual stars 
based, in turn, on stellar-evolution models or measurements in 
stellar binaries, assuming Newtonian dynamics. 

^ When we speak of SI, we are interested in what happens when 
we scale all lengths and times by a factor A. However, a the¬ 
ory usually involves other degrees of freedom, such as the gravi¬ 
tational potential, that are not directly measurable, call them 
Ip collectively. To these we can assign arbitrary scaling di¬ 
mensions, such that it is decreed that under the scaling of 
eq.(l), xp{r,t) —>■ A^-i/ifr/A, f/A). By multiplying ip by powers 
of the constants, we can always normalize them so that their 
[Z][i][m] dimensions match their scaling dimension. Namely, that 
if [ip] = [ifitr [m]**, then its ^ = /3 + 7 . With this choice of 
dimensions, simultaneous change of the length and time units 
by the same factor, under which all equations are invariant, is 
equivalent to a scaling transformation (which acts only on the 
degrees of freedom) followed by the change in the values of the 
constants. Thus SI is equivalent to invariance of the constants 
under the unit change. See more details in Ref. [12]. 
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general relativity)-summarized by a = — V(/), Acjj = 
AiiGp {(j) has zero scaling dimension), or by a ~ MGjr^- 
is not SI, is manifest by the appearance of G, whose value 
does change under the said change of units.® These the¬ 
ories can also not be modified into a SI theory without 
replacing G with other constants that are invariant to 
the change of units. 

In principle, several constants with different dimen¬ 
sions /3, 7 may appear in such a theory. In this first de¬ 
cision fork we assume, for parsimony, that only one such 
dimensioned constant appears. Such parsimony is not 
a mere whim. Reason advocates, and historical prece¬ 
dence with relativity and quantum theory concurs, that 
a major breakdown in well established physics does not 
introduce more than one new physical constant. 

Thus, G is banished from the DML, and replaced with 
another “gravitational constant”, call it Qo (dimension¬ 
less parameters may appear, in principle). 

We saw that SI predicts the asymptotic constancy of 
the rotational speed around an isolated body. Now that 
we accept that only one constant appears in the DML 
theory, the dependence of on the body’s mass is de¬ 
termined as follows: Since can depend only on the 
central mass, we must have, on dimensional grounds, 

oc (2) 

For /3 = 0, namely, when Qq has dimensions of velocity 
to some power, we get that does not depend on the 
mass and is a universal constant. This can be rejected at 
the outset because it clearly conflicts with observations. 

Since (3 ^ 0 we can assume, without loss of generality, 
that j3 = —1, which will be our choice of convenience 
henceforth.® With this standardized choice of units for 
Qo we get oc . It remains to fix the numer¬ 

ical normalization of Qo- We take it so that 

Ko = (QoM)i/T (3) 

(In MOND, Qo = ^o-) This mass-asymptotic-speed re¬ 
lation, M oc V2 is an important prediction of SI plus the 
assumption of only one new constant. 

The solitary appearance of Qo also leads to relations 
between different DML systems: An immediate corol¬ 
lary is that every NR, DML system is a member of a 
two-parameter family of homologous systems of different 
masses and sizes, with a unique relation between masses 
and velocities, independent of the size. 

To see this, start with any NR, self-gravitating, DML 
system, S. Make a change of units of mass, length, and 


® In standard units. Equivalently, we can work with units in which 
G = 1 , but then masses are not invariant. 

^ Starting with some Qo that has some 7 , and any ^ ^ 0, we 
could choose to work, instead, with Qo = Qq DML 

gravitational constant. This would have indeed set ^ ^ — 1 

for Qo (and change 7 to 7 = — 7 /^). 


time such that the values of attributes with these di¬ 
mensions change as: m —>■ = 6m^ i ^ 

and t ^ t' = The value of Qo remains the 

same in the new units, while the values of velocities be¬ 
come V —>■ v' = Since all equations are invariant 

to a change of units, a system homologous to S with 
the tagged attributes is also a solution of the DML the¬ 
ory (with the same value of Qo). We thus generate a 
whole two-parameter family of homologous solutions of 
all masses and sizes, within which all velocities (rotations, 
dispersions, etc.) are independent of size, and scale as 

y oc (QoM)i/T (4) 

The dimensionless proportionality factor in eq.(4) may 
be different for different families, and may depend on 
dimensionless attributes, such as the different mass ratios 
of constituents, velocity anisotropies, shape parameters, 
etc.^® 

For example, it follows that all DML-governed thin- 
disc galaxies with an exponential surface density profile, 
of all masses, and all exponential scale lengths, have the 
same value of /M. 

Despite the similar appearance, relation (4) is very dif¬ 
ferent from the M — VQ relation (3). The one concerns 
asymptotics, the other bulk properties; the one applies 
(asymptotically) to all systems, the other only to DML 
systems; the one predicts a universal equality with no 
scatter, the other not, as the ratio /M may depend 
on dimensionless attributes of the system. 

IV. THE UMBRELLA THEORY 

The DML theory, epitomized by the solitary appear¬ 
ance of Qo, is only part of the weak-field physics needed 
to describe gravitating systems. We know that it can 
only be a limit of a more complete theory, another limit 
of which is the weak-field general relativity (or Newto¬ 
nian dynamics), epitomized by the solitary appearance 
of G. The umbrella theory involves both constants. 

A. The boundary constant 

We expect the boundary between the validity domains 
of these two limits to be defined by some dimensioned 
constant, as h defines the boundary between the clas¬ 
sical and quantum domains, or c the classical-relativity 
boundary. Again, being loath to add yet more constants, 
I assume that this “boundary constant” is constructed 
from G and Qq. 

We want to retain the weak equivalence principle, 
namely the universality of free fall of test particles, in 


A different, two-parameter scaling, applies, of course, in Newto¬ 
nian dynamics due to the solitary appearance of G. 
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the DML, as well as in the full umbrella theory. This 
dictates that the “boundary constant” has no mass di¬ 
mensions. Otherwise, the question of whether some test 
particle motions are described by the DML or by New¬ 
tonian dynamics will depend on the mass of the particle, 
not only on the motion. 

To see this, consider a test particle of mass m in the 
system, and let B{tp, G, Qo) be a dimensionless criterion 
number, whose value tells us whether we are in the DML 
or not (similar to hu/kT in connection with the black- 
body spectrum in the classical-quantum case, or v/c and 
MG/rc^ in relativity). Here, f/' represents some local and 
nonlocal characteristic(s) of the particle’s orbit. Univer¬ 
sality of free fall dictates that il) cannot depend on m. 
Since B is dimensionless, it is invariant under a change 
of the mass units, under which the values of G and Qo 
do change. So, only the ratio qo = Qo/G, whose value 
does not change (since with our standardized choice of 
Qo it has the same mass dimensions as G), can appear in 
B. Furthermore, B{'ip, qo) can depend only on that are 
either dimensionless or of the dimensions of qo and then 
appear as ijj/qo- 

This is the case with c as the classical-relativistic 
“boundary constant”, which does not break universal¬ 
ity of free fall. On the other hand, h, as the classical- 
quantum “boundary constant”, does have mass dimen¬ 
sions, leading to manifest breakdown of the universality 
of free fall; For a given motion (velocity), the de Broglie 
wavelength of a particle depends on its mass.^^ The ap¬ 
pearance of G in Newtonian dynamics is not an obstacle, 
as it is not a boundary constant, only a conversion con¬ 
stant between inertial and gravitational masses. 

The boundary can thus be defined by qo, but it may 
depend on dimensionless characteristics of the orbit or 
phenomenon in question, such as the type of orbit (e.g., 
whether it is circular or not, etc.). So, for certain phe¬ 
nomena the actual boundary may, in principle, be a large 
or small multiple of go- I assume that this does not hap¬ 
pen, namely, the constant that emerges from the DML’s 
value of Qo (as qo = Qo/G) also marks the approxi¬ 
mate boundary for all phenomena. Such is the case, of 
course, for the transitions from classical physics to quan¬ 
tum physics and to relativity. The last verdict on this 
assumption is left to the observations. 

If all the test particles from which a system is made are 
always on one side of the boundary, we may say that the 
whole system is on that side. In this case we may express 


Some observational evidence for this comes, for example, from 
the fact that different test bodies, such as gas atoms, stars, etc. 
show the same rotation curves within observational errors. The 
bounds on departure from the principle are much looser than in 
laboratory and solar system measurements. 

So, for example, the analogue of the Bohr radius of a test mass 
m in the gravitational field of a finite mass M, tb ~ lm?MG, 
depends on m. (However, relative to the Planck length, Vp, we 
have tbItp M^/m?M, where Mp is the Planck mass.) 


the characteristic value of ip for a typical test particle in 
terms of the global properties of the system (such as mass 
and size), and get a global criterion for the system to be 
on that side of the boundary. For example, in the case of 
MOND where qo = Oq, such a system parameter may be 
(P jR, or MGjR?, or {MAoi^l'^jR. This is in the basis 
of the statement that a low-surface-density system: with 
MG/RP^ <C Oq is wholly in the DML. Note that here the 
total mass M does appear; this is not in conflict with the 
weak equivalence principle. 

The dimensions of qo are 

[go] = [£]'^-3[t]-(^-2). (5) 

In terms of go, eq. (3) reads 

= (goGM)i/^ (g) 

and eq.(4) reads 

U oc {qoGMf/\ (7) 

For example, if dynamics is modified beyond a certain 

scale length, lo, we have to take 7 = 2 , so go ~ I'q ^ 
is an inverse length. This predicts IQ oc If we 

want to modify dynamics below a critical frequency, wq, 
we have to take 7 = 8 , so go ~ Wg. In this case oc 
The MOND choice is 7 = 4, namely, go = Oq is an 
acceleration, which predicts oc 


B. A counterexample 

To recapitulate an important point: I have assumed 
that if the NR theory that I treat all along is attained as 
the NR limit of some relativistic theory, then this limit 
is strictly c —>■ 00 , so that c does not appear beside G in 
the NR limit. The NR theory is then assumed to have 
both a Newtonian limit, controlled, as usual, by G only, 
and a SI limit controlled by only one constant Qo- 

As an instructive exception I explain why so-called 
F{TZ) theories, or rather their NR limits, are not of the 
type discussed here, which disfavors them as candidate 
alternatives to dark matter. 

In such theories, the Einstein-Hilbert Lagrangian L oc 
AG~^TZ {TZ is the Ricci scalar), is replaced by L oc 
AG~^iQ‘^F{£QTZ), where £0 is some length constant of 
the theory. The strict NR limit, which assumes expan¬ 
sion near a Minkowski space time, is attained by taking 
c 00 , with £o/c = To and c^/G fixed. In this limit, 
£oTZ oc Tq A(/), where </> is the gravitational potential de¬ 
fined by goo = —I — 2(()/c^ (where goo is the time-time 
component of the metric). The Lagrangian then goes to 

Ljvr oc c^G“Vo“^F(to A(/)). ( 8 ) 

But this Ljvr, which involves both G and c, does not have 
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a Newtonian limit for any choice of F{x)-^^ (The particle 
Lagrangian is assumed intact and does not introduce any 
constants except masses.) 

The NR Lagrangian, does have a SI limit, which 
necessitates that in this limit F{x) oc so L^r oc 

with Qo = Gjclo, the dimensions of which 
are [Qo] = so we get a theory with 

7 = 1. This theory is conformally invariant in the three- 
dimensional Euclidean space. 

This higher-order theory (which suffers from the usual 
maladies of such theories) demonstrates also the need for 
our assumption that the asymptotic dynamics depends 
only on the total mass. The field equation 

A[(A())) 1 / 2 ] oc Qop, ( 9 ) 

has four independent spherical vacuum solutions: (j) = 
{QoMf X [(r/ro)"\ ln{r/ro), const., {r/rof]. Only 
the 2nd gives asymptotically flat RCs. That is because 
for the others solutions, the asymptotic dynamics, which 
is determined by V()), depends not only on the total mass, 
but also on structural attributes with the dimensions of 
length, which enter through rg, and which must also be 
scaled. So, the scaled orbit does not correspond to the 
same mass system. 


C. The discrepancy: definition and predictions 

The discrepancy between the observations and the pre¬ 
dictions of standard dynamics without dark matter is 
usually found and characterized by comparing the mea¬ 
sured acceleration, g, of test particles on some orbits at 
some location in a gravitating system, with the accel¬ 
eration, g^,, predicted by standard dynamics assuming 
that only the observed matter is present. Thus, as some 
measure of the discrepancy we can take 7 = g/g^, where 
9 = Is!) 9 n = IgArj- The ratio ry is sometimes called “the 
mass discrepancy”. The heuristic reason for this is clear: 
For a spherical system, Newtonian dynamics say that at 
a given radius the acceleration is proportional to the en¬ 
closed mass. So rj would indeed be the ratio between 
the enclosed mass needed to produce g to that produc¬ 
ing gjv. But in general, when the mass distributions of 
baryons and of the dark matter-real or fictitious-are not 
spherical, and are different, the acceleration does not gen¬ 
erally measure some enclosed mass, so the “acceleration 
discrepancy” ry is not really a “mass discrepancy”; the 
term can only be used loosely. In fact, since in general g 
and gjv are not parallel, we could define several accelera¬ 
tion ratios, for example using separately the acceleration 
component ratios. 


In this connection, note the following cautionary exam¬ 
ple: Ref. [13] analyzed the out-of-the-galactic-plane dy¬ 
namics in the Milky Way. As regards MOND, they con¬ 
sidered, separately, the ratios of the radial and perpendic¬ 
ular components of g and g^, call them ijr = 9 ^/ 9 ^^, and 
Vz = 9 ^/ 9 N, respectively. They plot these discrepancy 
measures against the acceleration. They find that, in ac¬ 
cordance with the predictions of MOND, both discrepan¬ 
cies set in roughly below an acceleration of 10 “®cm s“^, 
and that ry^, rjz increase with the decrease in g as MOND 
prescribes. There is no known reason why these clear 
systematics should happen in the dark-matter paradigm. 
But Ref. [13] fault MOND since they find that ry^. and 
r]z are not quite the same. The MOND expression they 
used erroneously, gjv = gg{g/ao), applies only in high- 
symmetry situation, such as spherical symmetry, or the 
axisymmetric case relevant for rotation curves. In such 
cases, clearly, g and g^ must be parallel from symmetry, 
and so automatically ry^ = Vz- However, there is no full- 
fledged MOND theory that predicts them to be the same 
for the analysis of Ref. [13]. It is also to be noted that the 
results of such an analysis depend on the adopted, but 
rather uncertain, model for the distribution of baryons in 
the Milky Way, which determines the values of gN(R). 

Furthermore, in Newtonian dynamics, the accelera¬ 
tions of all test particles are the same at a given posi¬ 
tion; so ry would be a function of position only. This is 
also the case in modified-gravity theories where the NR 
accelerations are still given by the gradient of a poten¬ 
tial. But in a subclass of modified-dynamics theories, so 
termed “modified inertia” theories, defined and described 
in Ref. [14], the accelerations may depend not only on 
location but on the whole orbit. So particles on different 
orbits might give rise to different values of ry even if their 
acceleration is measured at the same position. 

Returning to our general discussion, what do our 
modified-dynamics theories predict for ry? In general, we 
expect that ry should be correlated with the value of sys¬ 
tem attributes, q, that have the dimensions of qo, such 
that when q ^ q^,, g k, 1 (no discrepancy), while for 
g < Oo, ry > 1 .^"^ 

When g is determined from the asymptotic parts of 
rotation curves of disc galaxies, our theories predict 
g = V^lr, with VQ given by eq.( 6 ), while g^ = MGjr^, 
where M is the total (baryonic) mass of the galaxy. Thus, 


V = Vn 


9o 

9-R.{r)' 


( 10 ) 


where (/^(r) = V2~"^lr (subscript TZ indicates that the 
quantity pertains to rotation curve). q^is a local param¬ 
eter that has the dimensions of qo, according to eq. (5). 


In the relativistic theory, the general-relativity limit is obtained 
by taking F{x) ^ X? if we take this limit in the NR La¬ 
grangian, it becomes a complete derivative, and is unacceptable. 
The two limits—the NR limit and the GR/Newtonian limit—are 
not interchangeable. 


On the asymptotically flat section of rotation curves, the discrep¬ 
ancy increases with radius. Also, a quantity with the dimensions 
of go has dimensions of velocity/length; so must decrease as 
the inverse radius. Thus, this is indeed the sign of the inequali¬ 
ties. 
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The same result is gotten when 77 is obtained from weak 
lensing by galaxies, when analyzed consistently assuming 
asymptotically flat rotation curves [ 10 ]. 

More generally, defining qfi{r) = jr for the 

fnll rotation curve, we expect 77 ,^ to be strongly corre¬ 
lated with, if not quite a function of, x = q-,i(r)lqo. Any 
such correlation is predicted to give 77 « 1 for x ^ 1 , and 
rj = p/x ioT X where p ~ 1 is a proportionality fac¬ 
tor that may depend on the exact phenomenon and orbit, 
and other aspects entering the measurement of 77 . The 
latter proportionality stems from SI: From their defini¬ 
tion, 77 scales as while x scales as jr. Since V is 

invariant, we have under scaling 77 —> Ar;, while x —>■ x/A. 

When X ^ 1 is achieved on the asymptotic parts of 
rotation curve we have p = 1 due to our normalization 
of Qo and qo- But for other instances of x ^ 1 only 
proportionality is predicted. For example, circular, but 
non-asymptotic orbits, or radial orbits, or the use of dif¬ 
ferent acceleration components to define 77 , might give 
other values of p. These details depend on the theory. 

Modified inertia theories of the type discussed in Ref. 
[14] do predict a strict functional dependence^^ 

Vk = liq^/qo), /(a;>l)«il, /(x < 1 ) x"\ 

( 11 ) 

That /(x <C 1) « X ^ follows from SI and the standard¬ 
ized definition of qo, since, for the same reasons as above, 
77 scales like A, while q (hence x) scale as A“^. This 
functional dependence does not extend beyond rotation 
curves; in fact we do not have a general result for such 
“modified-inertia” theories. 

Now consider how 77 varies not within a given system, 
but across different systems. Write some global discrep¬ 
ancy measure, rjg, in terms of the global attributes of 
DML systems: the total mass, M, some measure of the 
radius, R (e.g., the scale length of exponential discs, the 
half-mass radius, etc.), and some representative measure 
of the internal speeds, a (e.g. the mass-weighted root- 
mean-squared velocity). 

r]g = {a'^/R)/{MG/R^) = a'^R/MG. (12) 

Also, define a global attribute with the dimensions of qo'. 

Qo = a-^-VR. (13) 

The scaling arguments around eq.(4) tell us that within 

each two-parameter family of systems, 

Vg=Kg — , (14) 

qg 

where Kg is a dimensionless constant that might depend 
on dimensionless attributes of the family, and on the ex¬ 
act definition of our global system attributes (unlike the 


Ref. [14] deals specifically with MOND, but the same argumen¬ 
tation applies to our more general 7 case. 


case of rotation curves where the choice of these param¬ 
eters is well defined). 

For Newtonian systems, which have qg ^ qo, the New¬ 
tonian virial relation implies r/g ~ 1 . 

Now, to compare between different families of DML 
systems, we use the assumption of no dimensionless con¬ 
stants that much differ from unity. This tells us that 
in each family there are members near the MOND- 
Newtonian boundary where we have both qg ^ qo and 
77 g ~ 1. This means that Kg ~ 1 for all the families. 

The important consequence is thus that pg is strongly 
correlated with qg for all systems, such that pg ~ qo/qg 
when qg qo, and pg ^ 1 when qg ^ qo- In other words, 
a strong correlation of the form described in eq.(ll) is 
predicted. 

For example, if the modification hinges on distance 
(7 = 2 ), introducing as boundary a new length ig, we 
predict a correlation^® 

P^f{R/io). (15) 

In the case of MOND (7 = 4 ), we predict a correlation 

V^fig/ao), (16) 

which can also be written as [ 1 ] 

V ~ v{9N/ao), (17) 

because 77 = glgN- For v{y) we have iy{y ^ 1) w 1, and 
in the DML, SI requirs u{y <C 1) « ■ 


V. WHY AN ACCELERATION (WHY MOND)? 

To complete the road to MOND we have to pinpoint 
7 and hence the full dimensions of qo and Qo- 

We have identified above at least four regularities in 
which 7 appears and from which it can be determined, 
(i) The mass-asymptotic-speed relation eqs.(3)(6). (ii) 
The general virial relation between internal velocities 
and total mass, eqs.(4)(7). (iii) The identification of the 
“boundary constant” and its dimensions, (iv) The di¬ 
mensions of the system attribute with which the discrep¬ 
ancy correlates. 

The mass-asymptotic-speed relation is closely related 
to the observed Tully-Fisher relation, which started life 
as a correlation between galaxy luminosity in some pho¬ 
tometric band, and some measure of the rotational speed 
(e.g., the 21-cm line width). Many versions of such a rela¬ 
tion have been presented in the literature, using various 
measures of the rotation speed and various luminosity 


This and the ones below are predicted correlations not exact 
functional relations (hence the use of the ~ sign). Different phe¬ 
nomena, or different kinds of orbits may show somewhat different 
behaviors. 
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measures. As has been stressed by MOND from its ad¬ 
vent, and as is clear from our discussion above, SI mod¬ 
ified dynamics require strictly that the velocity measure 
used is the asymptotic constant speed. It is also required 
(as emphasized, e.g., in Ref. [15]) that the other variable 
be the total baryonic mass, not a luminosity, which at 
best measures the stellar mass, while many galaxies can 
have a substantial fraction of their mass in gas. As a 
result of these MOND admonitions we have seen more 
appropriate versions of the Tully-Fisher relation. For ex¬ 
ample, Ref. [16] found that using the asymptotic speed 
for a sample of high-mass galaxies gave a tighter corre¬ 
lation than using other velocity measures. For the first 
they found a slope of 0.24 ±0.07.^^ Reference [17] tested 
the predicted MOND baryonic Tully-Fisher relation for 
galaxies whose mass is dominated by gas (to minimize 
the effect of light-to-mass conversion for the stars). The 
result of this analysis, and of one including all types of 
galaxies, is shown in Fig. 3 of Ref. [4]. It was found that 
the M — V^ slope is very tightly constrained to be 7 ~ 4. 
Slopes of 7 = 2 (large distance modihcation) and 7 = 8 
(low frequency modification) are clearly excluded.^® 

Some of the weak-lensing results of Ref. [9], can also be 
viewed as constraining 7 : When using logarithmic poten¬ 
tials (equivalent to asymptotically flat rotation curves) 
to fit their lensing data (binned into several luminos¬ 
ity bins spanning 2.5 orders of magnitude). Ref. [9] 
plot the galaxy luminosity vs. the equivalent of V^. 
They analyzed separately red (elliptical) galaxies, and 
blue (disc) galaxies and find for the respective slopes 
7-1 = 0.24 ± 0.03 and 7 -^ = 0.23 ± 0.03. 

The general mass-velocity, virial relation for large- 
discrepancy systems, while well consistent with 7 = 4 
(e.g.. Refs. [18, 19]) does not provide as tight a con¬ 
straint as others due to the caveats mentioned in Sec. 
IV. 

To top the argument. Fig. 5 of Ref. [4] (produced by 
Stacy McGaugh)^®, shows the discrepancy, as derived 
from the rotation curves of 73 disc galaxies, at many 
radii, r. rj is first plotted as function of r, to see if a 
correlation of the kind (15) exists, which 7 = 2 predicts. 
No such correlation appears. 

Then, this figure shows a test the predictions of the 
case 7 = 4 (MOND): Equation (17) [and so also of (16)] 
is tested by plotting 77 ,^ against g^- It was found that, 
indeed, a strong correlation exists, which is consistent, in 


Ref. [16] still used a luminosity measure, but it was the iC-band 
luminosity, which is thought to be a good measure of the stellar 
mass. And the gas contribution to the mass is small in the sample 
studied. Their sample included a relatively small range of galaxy 
masses; reducing the leverage on the slope. 

Even at the time MOND was put forth, it was possible to re¬ 
ject large-distance modifications (7 = 2 ) on the basis of existing 
limits on the Tully-Fisher slope [2]. 

This is an up-to-date version superseding previous similar anal¬ 
yses, e.g. in Refs. [20] [21] for disc galaxies, and in Ref. [22] for 
pressure-supported systems. 


fact, with no intrinsic scatter. 

In such a figure, Oo appears in three different roles rel¬ 
evant to our discussion above: (i) Some of the points in 
the region 5 ^ ^ Uo (or g <C Oq) correspond to asymp¬ 
totic regions of the rotation curves. The fact that they 
give rj w (ao/ffw)^^^ reflects the validity of the prediction 
of eq. (6). (ii) Other points in this region correspond to 
radii within the bulk of galaxies that show large mass dis¬ 
crepancies everywhere (so called “low-surface-brightness 
galaxies”). Here, the more general prediction is vindi¬ 
cated. (iii) The transition from ry = 1 to the DML asymp¬ 
totic regime occurs at ~ Oq of the same value. The tran¬ 
sition occurs in a span of accelerations of order Uq itself 
(roughly between ao/2 and 2ao). This vindicates our sur¬ 
mise that no new dimensionless constants are introduced, 
so that Oo defined from the DML is also the “boundary 
constant”, and the width of the transition region is of the 
order of Oq itself. 


VI. SUMMARY AND DISCUSSION 

I have considered a “purist” picture of modihed dy¬ 
namics as an alternative to dark matter. In this, only 
one new dimensioned constant is introduced by the new 
dynamics, which might be taken as Oq, or equivalently as 
Aq. I also assumed that no dimensionless constants that 
differ much from unity enter the new dynamics. Such 
was the case in the classical-to-quantum and classical-to- 
relativistic extensions. 

Here again are the main signposts on the road to 
MOND: 

(i) Asymptotic flatness of rotation curves-elevated to 
a status of an axiom-implies that within a given galaxy 
orbital times are proportional to the radius. Among dif¬ 
ferent large-discrepancy objects, velocities depend only 
on mass not on size; namely, dynamical times are propor¬ 
tional to size. I concluded that it behooves us to adopt 
SI dynamics for systems showing large mass discrepancies 
(so termed DML dynamics). 

(ii) I assumed that, replacing G, a single dimensioned 
gravitational constant, Qo, appears in the DML (apart 
from particle masses), which can be generally standard¬ 
ized to have dimensions [Qo] = [’7r]“^[^]^[t]“’^. 

(iii) I assumed that the Umbrella theory encompassing 
the DML and Newtonian dynamics, which involves G and 
Qo, does not involve additional dimensioned constants, 
nor dimensionless constants that much differ from unity. 

(iv) I assumed that the Umbrella theory satisfies the 
universality of free fall, and hence the constant that 
marks the boundary between the two limits of the theory 
has to be (after some standardization) Qq = Qq/G, which 
has dimensions [qo] = [£]^“®[t]“(^“^^. 

(v) Finally, based on many pieces of evidence I argued 
that 7 = 4, pinpointing Qo = Ao as the MOND gravita¬ 
tional constant, and Qq = Oq as the MOND acceleration 



constant. 

MOND, with its revolving around an acceleration, also 
turns out to have an added appeal in that the measured 
value of Oo might have strong cosmological connotations 
(e.g., Ref. [1]). For example, 27rao « cHq « c^(A/3)^/^, 
where Hq is the Hubble constant, and A the “cosmologi¬ 
cal constant”. There would be no such a connection had 
the transition from standard dynamics to the SI dynam¬ 
ics occurred at some critical lengths or time. 

The above coincidence also means that the MOND 
length, = c^/oo, is of the order of the Hubble distance, 
or the de Sitter radius associated with the cosmological 
constant. Thus, the role of Oo as a “boundary constant” 
can also be expressed as the role of the cosmologically- 
significant length as a boundary constant: There is a 
marked transition from non SI dynamics when the length 
(.{a) = (? ja, associated with an acceleration a, is smaller 
than Im: to SI dynamics when l'> Im [4, 23]. 

While my minimalist and purist starting point may, 
in the end, prove too restrictive in the case of MOND, 
I see no cogent reasons, at present, to relax it. On the 
other hand, there are quite a few theories that purport to 
underlie MOND and that transgress these axioms, such 
as hybrid, dark-matter-MOND theories, or those that al¬ 
low Oo to be system dependent (see, e.g.. Ref. [4] for a 
review). 

Note also, in this connection, that the salient aspects 
of MOND can be described as a configuration-dependent 
renormalization of Newton’s constant. For DML config¬ 
urations the main MOND results amount to a renormal¬ 
ization 

G ^ A/a = (ao/a)G = [^{a) /£m]G. (18) 

A mass discrepancy also appears in the context of the 
measured expansion history of The Universe after the 
decoupling of baryons from radiation, a period when 
NR matter strongly dominates over radiation. As evi¬ 
denced by several observations, within standard dynam¬ 
ics, baryons alone cannot account for the expansion his¬ 
tory. Results are standardly expressed in terms of the 
baryon and cosmological dark-matter densities. But we 
can say, alternatively, using the best figures we now have 
[24], that baryons fall short by a factor of 27r (to within 
the uncertainties, which are a few percents). Instead of 
interpreting this as a ^ 1 : 5 baryon-to-dark-matter ra¬ 
tio, in the spirit of MOND we would note, instead, that 
a renormalization G —>■« 27rG would have an equivalent 
effect.^^’^^ 


Of course we cannot say that the data gives strictly 7 = 4, since 
it is subject to measurement and systematic uncertainties. But 
certainly 7 = 4 is by far the most acceptable integer. 

Note that “dark matter” is also required in standard dynamics 
to explain structure formation and the microwave background 
power spectrum. 

For much earlier times, Ref. [27] derived a rather tight limit 
on the value G' of G that enters the expansion history at the 


The idea of configuration- or experiment-dependent 
renormalization of G in the basis of MOND is appeal¬ 
ing. But, in default of a microscopic theory that leads 
to such varied results, the idea remains a rather empty 
truism. 

All such, more baroque, attempts at underlying theo¬ 
ries come mainly to remedy the shortcoming of MOND 
in accounting for the cosmological mass discrepancies. 
But we do not necessarily need to go beyond the purist 
picture to account for the cosmological anomalies. Rel¬ 
ativistic MOND theories that embody the purist tenets 
(in their NR limits) may still turn out to be rich enough 
in structure and possibilities to extend the MOND con¬ 
cept to cosmology. Of the presently known such classes 
of theories we may mention BIMOND (bimetric MOND) 
[25], and nonlocal metric theories [26]. 

There is an illuminating, formal analogy of our discus¬ 
sion in this paper with the phenomenon of propagation 
of gravity waves in non-viscous fluids, free of surface ten¬ 
sion. For a large pool of fluid of constant depth h, in a 
gravitational acceleration field g, the dispersion relation 
is 

= gk ■ tanh(fc/i), (19) 

where ui is the wave (angular) frequency, and k = jkj, 
where k is the wave vector. 

To underscore the analogy, write eq. (19) in terms of 
the following quantities: Instead of the constants g and 
h, use the combinations: Cs = {gh)^/"^ (a speed), and 
G = c^h = gh"^, and instead of the variables uj and k, 
use the variables = Cs/w, which is a length, and the 
acceleration a = c^k. We then have 

G/rg = a ■tanh{a/g). (20) 

This has the form of a basic MOND relation, with the 
constants g, h, G, and cf = gG analogous, respectively, 
to tto, Im, G, and Aq in MOND, and the variables a 
and G/rl the analogs of the acceleration calculated in 
MOND, and the Newtonian acceleration G/A (the grav¬ 
itating mass is 1 in this analogy). 

The Newtonian limit of MOND is analogous to the 
deep-pool limit, h ^ k~^ {Im 00 ), which can be 
affected by taking formally g 0 (uo ^ 0), with all 
other quantities fixed. Then dynamics involve only the 
constant G, which is not invariant to unit scaling. As 
in Newtonian dynamics the resulting relation G/r'^ « a 
is not SI. The DML corresponds to the opposite limit: 

^ 00 (oq —>■ 00 ), with Cs fixed, so G —>■ 0 (G —>■ 0). The 
dispersion relation then involves only c® and is SI, as it 
takes the form c^/r^ « 


time of nucleosynthesis, a few minutes after the big bang, when 
radiation (relativistic matter) was dominant. From the observed 
^He abundance they deduced an upper limit of \G'/G—l\ < 0.13. 
The analogy is even closer in the heuristic picture described in 
Ref. [23], where the wavelength 2'Kk~^ corresponds to the Unruh 
wavelength associated with the acceleration a. 
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This two-dimensional description of the system is 
clearly only an approximate, effective account of the full 
three-dimensional hydrodynamical behavior. It is only 
valid in the limit where the bulk speed of sound of the 
fluid Cb ^ Cg. It not only exemplifies the MOND be¬ 
haviors in the two limits, but also demonstrates how an 
explicit “interpolating function” [/x(x) = tanh(a;)] be¬ 
tween the limits emerges from basis considerations in 
an effective theory (similar to what may happen in the 
case of MOND-as has been stressed many times). This 
interpolating function has, in fact, the exact asymp¬ 
totic behaviors, at both ends, as is required in MOND 
[/x(x » I) « I, ^{x ^ 1) « x\. The system also satisfies 
our other assumptions of a single dimensioned constant 
in each limit, and of no dimensionless constants that dif¬ 
fer much from unity. 

Such “mechanistic” analogies might inspire construc¬ 
tion of MOND theories (as in the case of the molecular- 
vortices auxiliary picture for the Maxwell equations). 

And, add a footnote with a more formal explanation: 
For a test particle of mass m, let I?('0, G, Qo) be a dimen¬ 
sionless “criterion number”, whose value tells us whether 
we are in the DML or not (similar to hv/kT in connection 

Note that here the total mass M does appear; this is not in 

conflict with the weak equivalence principle. 


with the BB classical-quantum case, or v/c in relativity). 
Here, ip represents some characteristic(s) of the particle. 
Universality of free fall dictates that ip cannot depend on 
TO. Then, since B is dimensionless, it is invariant under a 
change of the mass units, under which the values of G and 
Qo do change. So, only the ratio Qo = Qo/G, whose value 
does not change, can appear in B. Furthermore, B{ip,qo) 
can depend only on ip that are either dimensionless or of 
the dimensions of qo and appear as ip/qo- The boundary 
is thus defined by qo, but may depend on dimensionless 
characteristics, such as the type of orbit. 

If all the test particles from which a system is made 
are always on one side of the boundary we may say that 
the whole system is on that side. In this case we may ex¬ 
press the characteristic value of ip for a characteristic test 
particle in terms of the global properties of the system 
(such as mass and size), and get a global criterion for the 
system to be on that side of the boundary. For example, 
in the case of MOND where qo = Uq, such a system pa¬ 
rameter may be cP' jR, or MGjRP, or {MAqYPIR. This 
is in the basis of the statement that a low-surface-density 
system: with MG/R^ <^ao is wholly in the DML.^^ 
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